Systematic in vitro analysis of therapy resistance in glioblastoma cell lines by integration of clonogenic survival data with multi-level molecular data.

Despite intensive basic scientific, translational, and clinical efforts in the last decades, glioblastoma remains a devastating disease with a highly dismal prognosis. Apart from the implementation of temozolomide into the clinical routine, novel treatment approaches have largely failed, emphasizing the need for systematic examination of glioblastoma therapy resistance in order to identify major drivers and thus, potential vulnerabilities for therapeutic intervention. Recently, we provided proof-of-concept for the systematic identification of combined modality radiochemotherapy treatment vulnerabilities via integration of clonogenic survival data upon radio(chemo)therapy with low-density transcriptomic profiling data in a panel of established human glioblastoma cell lines. Here, we expand this approach to multiple molecular levels, including genomic copy number, spectral karyotyping, DNA methylation, and transcriptome data. Correlation of transcriptome data with inherent therapy resistance on the single gene level yielded several candidates that were so far underappreciated in this context and for which clinically approved drugs are readily available, such as the androgen receptor (AR). Gene set enrichment analyses confirmed these results, and identified additional gene sets, including reactive oxygen species detoxification, mammalian target of rapamycin complex 1 (MTORC1) signaling, and ferroptosis/autophagy-related regulatory circuits to be associated with inherent therapy resistance in glioblastoma cells. To identify pharmacologically accessible genes within those gene sets, leading edge analyses were performed yielding candidates with functions in thioredoxin/peroxiredoxin metabolism, glutathione synthesis, chaperoning of proteins, prolyl hydroxylation, proteasome function, and DNA synthesis/repair. Our study thus confirms previously nominated targets for mechanism-based multi-modal glioblastoma therapy, provides proof-of-concept for this workflow of multi-level data integration, and identifies novel candidates for which pharmacological inhibitors are readily available and whose targeting in combination with radio(chemo)therapy deserves further examination. In addition, our study also reveals that the presented workflow requires mRNA expression data, rather than genomic copy number or DNA methylation data, since no stringent correlation between these data levels could be observed. Finally, the data sets generated in the present study, including functional and multi-level molecular data of commonly used glioblastoma cell lines, represent a valuable toolbox for other researchers in the field of glioblastoma therapy resistance.

Transcriptomic landscape of radiation-induced murine thyroid proliferative lesions.

Thyroid carcinoma incidence rates in western societies are among the fastest rising, compared to all malignant tumors over the past two decades. While risk factors such as age and exposure to ionizing radiation are known, early-state carcinogenic processes or pre-lesions are poorly understood or unknown. This study aims at the identification and characterization of early-state radiation-associated neoplastic processes by histologic and transcriptomic analyses of thyroid tissues derived from a mouse model. Comprehensive histological examination of 246 thyroids (164 exposed, 82 non-exposed) was carried out. Proliferative and normal tissues from exposed cases and normal tissue from non-exposed cases were collected by laser-capture microdissection, followed by RNAseq transcriptomic profiling using a low input 3'-library preparation protocol, differential gene expression analysis and functional association by gene set enrichment analysis. Nine exposed samples exhibited proliferative lesions, while none of the non-exposed samples showed histological abnormalities, indicating an association of ionizing radiation exposure with histological abnormalities. Activated immune response signaling and deregulated metabolic processes were observed in irradiated tissue with normal histology compared to normal tissue from non-exposed samples. Proliferative lesions compared to corresponding normal tissues showed enrichment for mainly proliferation-associated gene sets. Consistently, proliferative lesion samples from exposed mice showed elevated proliferation-associated signaling and deregulated metabolic processes compared to normal samples from non-exposed mice. Our findings suggest that a molecular deregulation may be detectable in histologically normal thyroid tissues and in early proliferative lesions in the frame of multi-step progression from irradiated normal tissue to tumorous lesions.

Transcriptome network of the papillary thyroid carcinoma radiation marker CLIP2.

BACKGROUND: We present a functional gene association network of the CLIP2 gene, generated by de-novo reconstruction from transcriptomic microarray data. CLIP2 was previously identified as a potential marker for radiation induced papillary thyroid carcinoma (PTC) of young patients in the aftermath of the Chernobyl reactor accident. Considering the rising thyroid cancer incidence rates in western societies, potentially related to medical radiation exposure, the functional characterization of CLIP2 is of relevance and contributes to the knowledge about radiation-induced thyroid malignancies. METHODS: We generated a transcriptomic mRNA expression data set from a CLIP2-perturbed thyroid cancer cell line (TPC-1) with induced CLIP2 mRNA overexpression and siRNA knockdown, respectively, followed by gene-association network reconstruction using the partial correlation-based approach GeneNet. Furthermore, we investigated different approaches for prioritizing differentially expressed genes for network reconstruction and compared the resulting networks with existing functional interaction networks from the Reactome, Biogrid and STRING databases. The derived CLIP2 interaction partners were validated on transcript and protein level. RESULTS: The best reconstructed network with regard to selection parameters contained a set of 20 genes in the 1st neighborhood of CLIP2 and suggests involvement of CLIP2 in the biological processes DNA repair/maintenance, chromosomal instability, promotion of proliferation and metastasis. Peptidylprolyl Isomerase Like 3 (PPIL3), previously identified as a potential direct interaction partner of CLIP2, was confirmed in this study by co-expression at the transcript and protein level. CONCLUSION: In our study we present an optimized preselection approach for genes subjected to gene-association network reconstruction, which was applied to CLIP2 perturbation transcriptome data of a thyroid cancer cell culture model. Our data support the potential carcinogenic role of CLIP2 overexpression in radiation-induced PTC and further suggest potential interaction partners of the gene.

Comparative Transcriptomic Analysis of Temozolomide Resistant Primary GBM Stem-Like Cells and Recurrent GBM Identifies Up-Regulation of the Carbonic Anhydrase CA2 Gene as Resistance Factor.

About 95% of patients with Glioblastoma (GBM) show tumor relapse, leaving them with limited therapeutic options as recurrent tumors are most often resistant to the first line chemotherapy standard Temozolomide (TMZ). To identify molecular pathways involved in TMZ resistance, primary GBM Stem-like Cells (GSCs) were isolated, characterized, and selected for TMZ resistance in vitro. Subsequently, RNA sequencing analysis was performed and revealed a total of 49 differentially expressed genes (|log2-fold change| > 0.5 and adjusted p-value < 0.1) in TMZ resistant stem-like cells compared to their matched DMSO control cells. Among up-regulated genes, we identified carbonic anhydrase 2 (CA2) as a candidate gene correlated with glioma malignancy and patient survival. Notably, we describe consistent up-regulation of CA2 not only in TMZ resistant GSCs on mRNA and protein level, but also in patient-matched clinical samples of first manifest and recurrent tumors. Co-treatment with the carbonic anhydrase inhibitor Acetazolamid (ACZ) sensitized cells to TMZ induced cell death. Cumulatively, our findings illustrate the potential of CA2 as a chemosensitizing target in recurrent GBM and provide a rationale for a therapy associated inhibition of CA2 to overcome TMZ induced chemoresistance.

A prognostic mRNA expression signature of four 16q24.3 genes in radio(chemo)therapy-treated head and neck squamous cell carcinoma (HNSCC).

Previously, we have shown that copy number gain of the chromosomal band 16q24.3 is associated with impaired clinical outcome of radiotherapy-treated head and neck squamous cell carcinoma (HNSCC) patients. We set out to identify a prognostic mRNA signature from genes located on 16q24.3 in radio(chemo)therapy-treated HNSCC patients of the TCGA (The Cancer Genome Atlas, n = 99) cohort. We applied stepwise forward selection using expression data of 41 16q24.3 genes. The resulting optimal Cox-proportional hazards regression model included the genes APRT, CENPBD1, CHMP1A, and GALNS. Afterward, the prognostic value of the classifier was confirmed in an independent cohort of HNSCC patients treated by adjuvant radio(chemo)therapy (LMU-KKG cohort). The signature significantly differentiated high- and low-risk patients with regard to overall survival (HR = 2.01, 95% CI 1.10-3.70; P = 0.02125), recurrence-free survival (HR = 1.84, 95% CI 1.01-3.34; P = 0.04206), and locoregional recurrence-free survival (HR = 1.87, 95% CI 1.03-3.40; P = 0.03641). The functional impact of the four signature genes was investigated after reconstruction of a gene association network from transcriptome data of the TCGA HNSCC cohort using a partial correlation approach. Subsequent pathway enrichment analysis of the network neighborhood (first and second) of the signature genes suggests involvement of HNSCC-associated signaling pathways such as apoptosis, cell cycle, cell adhesion, EGFR, JAK-STAT, and mTOR. Furthermore, a detailed analysis of the first neighborhood revealed a cluster of co-expressed genes located on chromosome 16q, substantiating the impact of 16q24.3 alterations in poor clinical outcome of HNSCC. The reported gene expression signature represents a prognostic marker in HNSCC patients following postoperative radio(chemo)therapy.

A genomic copy number signature predicts radiation exposure in post-Chernobyl breast cancer.

Breast cancer is the second leading cause of cancer death among women worldwide and besides life style, age and genetic risk factors, exposure to ionizing radiation is known to increase the risk for breast cancer. Further, DNA copy number alterations (CNAs), which can result from radiation-induced double-strand breaks, are frequently occurring in breast cancer cells. We set out to identify a signature of CNAs discriminating breast cancers from radiation-exposed and non-exposed female patients. We analyzed resected breast cancer tissues from 68 exposed female Chernobyl clean-up workers and evacuees and 68 matched non-exposed control patients for CNAs by array comparative genomic hybridization analysis (aCGH). Using a stepwise forward-backward selection approach a non-complex CNA signature, that is, less than ten features, was identified in the training data set, which could be subsequently validated in the validation data set (p value < 0.05). The signature consisted of nine copy number regions located on chromosomal bands 7q11.22-11.23, 7q21.3, 16q24.3, 17q21.31, 20p11.23-11.21, 1p21.1, 2q35, 2q35, 6p22.2. The signature was independent of any clinical characteristics of the patients. In all, we identified a CNA signature that has the potential to allow identification of radiation-associated breast cancer at the individual level.

Expression of miRNA-26b-5p and its target TRPS1 is associated with radiation exposure in post-Chernobyl breast cancer.

Ionizing radiation is a well-recognized risk factor for the development of breast cancer. However, it is unknown whether radiation-specific molecular oncogenic mechanisms exist. We investigated post-Chernobyl breast cancers from radiation-exposed female clean-up workers and nonexposed controls for molecular changes. Radiation-associated alterations identified in the discovery cohort (n = 38) were subsequently validated in a second cohort (n = 39). Increased expression of hsa-miR-26b-5p was associated with radiation exposure in both of the cohorts. Moreover, downregulation of the TRPS1 protein, which is a transcriptional target of hsa-miR-26b-5p, was associated with radiation exposure. As TRPS1 overexpression is common in sporadic breast cancer, its observed downregulation in radiation-associated breast cancer warrants clarification of the specific functional role of TRPS1 in the radiation context. For this purpose, the impact of TRPS1 on the transcriptome was characterized in two radiation-transformed breast cell culture models after siRNA-knockdown. Deregulated genes upon TRPS1 knockdown were associated with DNA-repair, cell cycle, mitosis, cell migration, angiogenesis and EMT pathways. Furthermore, we identified the interaction partners of TRPS1 from the transcriptomic correlation networks derived from gene expression data on radiation-transformed breast cell culture models and sporadic breast cancer tissues provided by the TCGA database. The genes correlating with TRPS1 in the radiation-transformed breast cell lines were primarily linked to DNA damage response and chromosome segregation, while the transcriptional interaction partners in the sporadic breast cancers were mostly associated with apoptosis. Thus, upregulation of hsa-miR-26b-5p and downregulation of TRPS1 in radiation-associated breast cancer tissue samples suggests these molecules representing radiation markers in breast cancer.

Genomic amplification of Fanconi anemia complementation group A (FancA) in head and neck squamous cell carcinoma (HNSCC): Cellular mechanisms of radioresistance and clinical relevance.

Radio (chemo) therapy is a crucial treatment modality for head and neck squamous cell carcinoma (HNSCC), but relapse is frequent, and the underlying mechanisms remain largely elusive. Therefore, novel biomarkers are urgently needed. Previously, we identified gains on 16q23-24 to be associated with amplification of the Fanconi anemia A (FancA) gene and to correlate with reduced progression-free survival after radiotherapy. Here, we analyzed the effects of FancA on radiation sensitivity in vitro, characterized the underlying mechanisms, and evaluated their clinical relevance. Silencing of FancA expression in HNSCC cell lines with genomic gains on 16q23-24 resulted in significantly impaired clonogenic survival upon irradiation. Conversely, overexpression of FancA in immortalized keratinocytes conferred increased survival accompanied by improved DNA repair, reduced accumulation of chromosomal translocations, but no hyperactivation of the FA/BRCA-pathway. Downregulation of interferon signaling as identified by microarray analyses, enforced irradiation-induced senescence, and elevated production of the senescence-associated secretory phenotype (SASP) appeared to be candidate mechanisms contributing to FancA-mediated radioresistance. Data of the TCGA HNSCC cohort confirmed the association of gains on 16q24.3 with FancA overexpression and impaired overall survival. Importantly, transcriptomic alterations similar to those observed upon FancA overexpression in vitro strengthened the clinical relevance. Overall, FancA amplification and overexpression appear to be crucial for radiotherapeutic failure in HNSCC.

Integration of a radiation biomarker into modeling of thyroid carcinogenesis and post-Chernobyl risk assessment.

Strong evidence for the statistical association between radiation exposure and disease has been produced for thyroid cancer by epidemiological studies after the Chernobyl accident. However, limitations of the epidemiological approach in order to explore health risks especially at low doses of radiation appear obvious. Statistical fluctuations due to small case numbers dominate the uncertainty of risk estimates. Molecular radiation markers have been searched extensively to separate radiation-induced cancer cases from sporadic cases. The overexpression of the CLIP2 gene is the most promising of these markers. It was found in the majority of papillary thyroid cancers (PTCs) from young patients included in the Chernobyl tissue bank. Motivated by the CLIP2 findings we propose a mechanistic model which describes PTC development as a sequence of rate-limiting events in two distinct paths of CLIP2-associated and multistage carcinogenesis. It integrates molecular measurements of the dichotomous CLIP2 marker from 141 patients into the epidemiological risk analysis for about 13 000 subjects from the Ukrainian-American cohort which were exposed below age 19 years and were put under enhanced medical surveillance since 1998. For the first time, a radiation risk has been estimated solely from marker measurements. Cross checking with epidemiological estimates and model validation suggests that CLIP2 is a marker of high precision. CLIP2 leaves an imprint in the epidemiological incidence data which is typical for a driver gene. With the mechanistic model, we explore the impact of radiation on the molecular landscape of PTC. The model constitutes a unique interface between molecular biology and radiation epidemiology.

Natural Cubic Spline Regression Modeling Followed by Dynamic Network Reconstruction for the Identification of Radiation-Sensitivity Gene Association Networks from Time-Course Transcriptome Data.

Gene expression time-course experiments allow to study the dynamics of transcriptomic changes in cells exposed to different stimuli. However, most approaches for the reconstruction of gene association networks (GANs) do not propose prior-selection approaches tailored to time-course transcriptome data. Here, we present a workflow for the identification of GANs from time-course data using prior selection of genes differentially expressed over time identified by natural cubic spline regression modeling (NCSRM). The workflow comprises three major steps: 1) the identification of differentially expressed genes from time-course expression data by employing NCSRM, 2) the use of regularized dynamic partial correlation as implemented in GeneNet to infer GANs from differentially expressed genes and 3) the identification and functional characterization of the key nodes in the reconstructed networks. The approach was applied on a time-resolved transcriptome data set of radiation-perturbed cell culture models of non-tumor cells with normal and increased radiation sensitivity. NCSRM detected significantly more genes than another commonly used method for time-course transcriptome analysis (BETR). While most genes detected with BETR were also detected with NCSRM the false-detection rate of NCSRM was low (3%). The GANs reconstructed from genes detected with NCSRM showed a better overlap with the interactome network Reactome compared to GANs derived from BETR detected genes. After exposure to 1 Gy the normal sensitive cells showed only sparse response compared to cells with increased sensitivity, which exhibited a strong response mainly of genes related to the senescence pathway. After exposure to 10 Gy the response of the normal sensitive cells was mainly associated with senescence and that of cells with increased sensitivity with apoptosis. We discuss these results in a clinical context and underline the impact of senescence-associated pathways in acute radiation response of normal cells. The workflow of this novel approach is implemented in the open-source Bioconductor R-package splineTimeR.

Transcriptomic analyses of the radiation response in head and neck squamous cell carcinoma subclones with different radiation sensitivity: time-course gene expression profiles and gene association networks.

BACKGROUND: Acquired and inherent radioresistance of tumor cells is related to tumor relapse and poor prognosis - not only in head and neck squamous cell carcinoma (HNSCC). The underlying molecular mechanisms are largely unknown. Therefore, systemic in-depth analyses are needed to identify key regulators of radioresistance. In the present study, subclones of the CAL-33 HNSCC cell line with different radiosensitivity were analyzed to identify signaling pathways related to the different phenotypes. METHODS: Subclones with altered radiosensitivity were generated by fractionated irradiation of the parental CAL-33 cells. Differences in radiosensitivity were confirmed in colony formation assays. Selected subclones were characterized at the genomic and transcriptomic level by SKY, array CGH, and mRNA-microarray analyses. Time-course gene expression analyses upon irradiation using a natural cubic spline regression model identified temporally differentially expressed genes. Moreover, early and late responding genes were identified. Gene association networks were reconstructed using partial correlation. The Reactome pathway database was employed to conduct pathway enrichment analyses. RESULTS: The characterization of two subclones with enhanced radiation resistance (RP) and enhanced radiosensitivity (SP) revealed distinct genomic and transcriptomic changes compared to the parental cells. Differentially expressed genes after irradiation shared by both subclones pointed to important pathways of the early and late radiation response, including senescence, apoptosis, DNA repair, Wnt, PI3K/AKT, and Rho GTPase signaling. The analysis of the most important nodes of the gene association networks revealed pathways specific to the radiation response in different phenotypes of radiosensitivity. Exemplarily, for the RP subclone the senescence-associated secretory phenotype (SASP) together with GPCR ligand binding were considered as crucial. Also, the expression of endogenous retrovirus ERV3-1in response to irradiation has been observed, and the related gene association networks have been identified. CONCLUSIONS: Our study presents comprehensive gene expression data of CAL-33 subclones with different radiation sensitivity. The resulting networks and pathways associated with the resistant phenotype are of special interest and include the SASP. The radiation-associated expression of ERV3-1 also appears highly attractive for further studies of the molecular mechanisms underlying acquired radioresistance. The identified pathways may represent key players of radioresistance, which could serve as potential targets for molecularly designed, therapeutical intervention.

Integrative analysis of the microRNA-mRNA response to radiochemotherapy in primary head and neck squamous cell carcinoma cells.

BACKGROUND: Head and neck squamous cell carcinoma (HNSCC) is a very heterogeneous disease resulting in huge differences in the treatment response. New individualized therapy strategies including molecular targeting might help to improve treatment success. In order to identify potential targets, we developed a HNSCC radiochemotherapy cell culture model of primary HNSCC cells derived from two different patients (HN1957 and HN2092) and applied an integrative microRNA (miRNA) and mRNA analysis in order to gain information on the biological networks and processes of the cellular therapy response. We further identified potential target genes of four therapy-responsive miRNAs detected previously in the circulation of HNSCC patients by pathway enrichment analysis. RESULTS: The two primary cell cultures differ in global copy number alterations and P53 mutational status, thus reflecting heterogeneity of HNSCC. However, they also share many copy number alterations and chromosomal rearrangements as well as deregulated therapy-responsive miRNAs and mRNAs. Accordingly, six common therapy-responsive pathways (direct P53 effectors, apoptotic execution phase, DNA damage/telomere stress induced senescence, cholesterol biosynthesis, unfolded protein response, dissolution of fibrin clot) were identified in both cell cultures based on deregulated mRNAs. However, inflammatory pathways represented an important part of the treatment response only in HN1957, pointing to differences in the treatment responses of the two primary cultures. Focused analysis of target genes of four therapy-responsive circulating miRNAs, identified in a previous study on HNSCC patients, revealed a major impact on the pathways direct P53 effectors, the E2F transcription factor network and pathways in cancer (mainly represented by the PTEN/AKT signaling pathway). CONCLUSIONS: The integrative analysis combining miRNA expression, mRNA expression and the related cellular pathways revealed that the majority of radiochemotherapy-responsive pathways in primary HNSCC cells are related to cell cycle, proliferation, cell death and stress response (including inflammation). Despite the heterogeneity of HNSCC, the two primary cell cultures exhibited strong similarities in the treatment response. The findings of our study suggest potential therapeutic targets in the E2F transcription factor network and the PTEN/AKT signaling pathway.

Genomic copy number analysis of Chernobyl papillary thyroid carcinoma in the Ukrainian-American Cohort.

One of the major consequences of the 1986 Chernobyl reactor accident was a dramatic increase in papillary thyroid carcinoma (PTC) incidence, predominantly in patients exposed to the radioiodine fallout at young age. The present study is the first on genomic copy number alterations (CNAs) of PTCs of the Ukrainian-American cohort (UkrAm) generated by array comparative genomic hybridization (aCGH). Unsupervised hierarchical clustering of CNA profiles revealed a significant enrichment of a subgroup of patients with female gender, long latency (>17 years) and negative lymph node status. Further, we identified single CNAs that were significantly associated with latency, gender, radiation dose and BRAF V600E mutation status. Multivariate analysis revealed no interactions but additive effects of parameters gender, latency and dose on CNAs. The previously identified radiation-associated gain of the chromosomal bands 7q11.22-11.23 was present in 29% of cases. Moreover, comparison of our radiation-associated PTC data set with the TCGA data set on sporadic PTCs revealed altered copy numbers of the tumor driver genes NF2 and CHEK2. Further, we integrated the CNA data with transcriptomic data that were available on a subset of the herein analyzed cohort and did not find statistically significant associations between the two molecular layers. However, applying hierarchical clustering on a 'BRAF-like/RAS-like' transcriptome signature split the cases into four groups, one of which containing all BRAF-positive cases validating the signature in an independent data set.

Dose-dependent expression of CLIP2 in post-Chernobyl papillary thyroid carcinomas.

A previous study on papillary thyroid carcinomas (PTC) in young patients who were exposed to (131)iodine from the Chernobyl fallout revealed an exclusive gain of chromosomal band 7q11.23 in exposed cases compared to an age-matched control cohort. CLIP2, a gene located within band 7q11.23 was shown to be differentially expressed between exposed and non-exposed cases at messenger RNA and protein level. Therefore, a standardized procedure for CLIP2 typing of PTCs has been developed in a follow-up study. Here we used CLIP2 typing data on 117 post-Chernobyl PTCs from two cohorts of exposed patients with individual dose estimates and 24 non-exposed controls to investigate a possible quantitative dose-response relationship of the CLIP2 marker. The 'Genrisk-T' cohort consisted of 45 PTCs and the 'UkrAm' cohort of 72 PTCs. Both cohorts differed in mean dose (0.59 Gy Genrisk-T, 1.2 Gy UkrAm) and mean age at exposure (AaE) (2 years Genrisk-T, 8 years UkrAm), whilst the median latency (16 years Genrisk-T, 18 years UkrAm) was comparable. We analyzed the association between the binary CLIP2 typing and continuous thyroid dose with logistic regression. A clear positive dose-response relationship was found for young PTC cases [age at operation (AaO) < 20 years, AaE < 5 years]. In the elder age group a higher proportion of sporadic tumors is assumed due to a negligible dose response, suggesting different molecular mechanisms in sporadic and radiation-induced cases. This is further supported by the association of elder patients (AaO > 20 years) with positivity for BRAF V600E mutation.